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Abstract. Intraday Variability of compact extragalactic radio sources can be interpreted as quenched scintillation 
due to turbulent density fluctuations of the nearby ionized interstellar medium. We demonstrate that the statistical 
analysis of IDV time series contains both information about sub-structure of the source on the scale of several 10 
fias and about the turbulent state of the ISM. The source structure and ISM properties cannot be disentangled 
using IDV observations alone. A comparison with the known morphology of the 'local bubble' and the turbulent 
ISM from pulsar observations constrains possible source models. We further argue that earth orbit synthesis fails 
for non- stationary relativistic sources and no reliable 2D-Fourier reconstruction is possible. 
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Intraday Variability (IDV) (see Wagner & Witzel ( pL995[ ) 
for a review) is a common phenomenon of flat-spectrum 
radio cores in quasars and BL Lacs. From light travel time 
arguments the observed brightness temperature of IDV 
sources are in the range of 10 16 - - 21 K and far in excess of 
the inverse Compton limit. The suggested intrinsic expla- 
nation require cither extreme Doppler boosting or special 
source geometries. Gravitational micro-lensing and scin- 
tillation in the ISM of our galaxy are also discussed as 
possible propagation effects causing variability. 

Two classes of IDV sources^ are distinguished accord- 
ing to their structure functions (<SF). One shows a contin- 
uous increase of SF with time-lag, while the other reaches 
a 'plateau' at a well defined time-scale tiuv (see Fig. 1). 

In this paper we will discuss the scintillation hypoth- 
esis for the 'plateau' class. Scintillation is caused by tur- 
bulence in the ISM, which is a stochastic process and the 
resulting time series are best analysed in terms of first 
order structure functions. Following Simonetti, Cordes & 
Heeschen (19851) we define the discrete structure function 



for time series f(ti) by 



SF( Tj ) = N t 
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.y is the normalisation; w(i) are weighting func- 



Here N i} 

tions, so that w(i)w(i + j) > 0, if a measurement at U 
and another measurement at t% + Tj were obtained. The 
weighting function also accounts for birring of unevenly 
sampled data. The structure function is related to the au- 
tocorrelation p(r) by SF(t) = 2[p(0) - p(r)]. 



1 Besides normal IDV two extre mely fast sources PKS 
0405-385 (Kedziora-Chudczer et al. |L997|) and J1819+3845 



(Dennett-Thorpe & de Bruyn 2000 ) have been found with very 
short time-scales of less than 1 hour. 
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Fig. 1. Structure function for IDV in total intensity for 
the quasar 0917+624. The SF oc t° has a slope of a w 
1, a 'plateau' with SF = 2m 2 for long time-lags and a 
characteristic timescale for decorrelation tiny at the first 
maximum. 



2. Quenched Scintillation 

Turbulence in the ionized ISM of our Galaxy produces 
density fluctuations of free electrons and the correspond- 
ing refractive index for radio wave propagation on a wide 
range of spatial scales (Armstrong, Rickett & Spangler 
1995). The power spectrum of density fluctuations in the 
ISM at wave number q along the line of sight (coordi- 
nate z) towards the quasar is described by a power law 
&(z,q) = Cf f (z)q~ t3 between a lower and upper cut-off 
<Ziow < q < Shigh- The lower cut-off in q corresponds to the 
driving wavelength of the turbulent cascade. The inertial 
range described by the power law ends at (?high, where 
energy is dissipated on small spatial scales. 

One distinguishes refractive and diffractive scattering 
at low frequencies and weak scattering above a ch arac- 
teristic transition frequency vt (e.g. Narayan 1992] ). For 
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quasars at high galactic latitude, t he tra nsition is expected 
between 1 and 5 GHz (e.g. Walker 1998). The angular scale 
for destructive interference of wavefronts due to optical 
path length differences is the Fresnel scale 



brightness V(q). The argument of the sinusoidal Fresnel 
term 



6 F = \ = 0.26mas 

Z7TZ 



GHz pc 



-1/2 



(2) 



Large fluctuations on scales smaller than Op leads to in- 
terference and diffractive scintillation, while fluctuations 
on scales larger than Op act as gradients in the refractive 
index and causes focusing and defocusing of light rays. 

Due to the 'large' size of incoherent extragalactic syn- 
chrotron sources 
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diffractive effects are unimportant, because the source size 
is always larger than the Fresnel scale for scattering in the 
ISM. Here z is the cosmological redshift, Ti, the brightness 
temperature, F v the observed flux density of the source, 
v the observing frequency, and S the Doppler boosting 
factor. 

For a scattering screen lying beyond the solar system 
and at frequencies v > 1 GHz the source size is always 
much larger than the Fresnel scale and the amplitude of 
variations is quenched. The amplitude of variations is re- 
duced and the typical variability time-scale is stretched by 
the large size of the source. 

The fluctuations in the ISM at scale q change with a 
characteristic turbulent velocity v(q), which is assumed 
much smaller than the relative velocity of earth and ISM. 
The scintillation pattern of intensity variations appears to 
be frozen in the ISM and the spatial variations are scanned 
by earth motion. The relative velocity of earth and ISM 
projected onto the sky has contributions from earth orbital 
motion, fEarth, the velocity of the solar system in the local 
standard of rest, wlsr.j and the peculiar velocities of clouds 
and sheets in the ISM, w c ioud, like those in the wall of the 
local bubble 



WEarth + «LSR + «cloud 



(4) 



Due to orbital motion, the velocity, v, traces an off-centred 
ellipse in the plane of the sky. The velocity transfers vari- 
ations on spatial scales x to temporal variations t — x/v 
and leads to an ann ual m odulation of the variability time- 
scale (Rickett et al. |200l|) . 



From theoretical considerations (Coles et al. 1987) the 
autocorrelation p(r) of light curves due to quenched scin- 
tillation can be calculated as an integral along the line of 
sight through the scattering medium of a Fourier trans- 
form q — > x 

p{x) = 87rr 2 e X 2 J dzFT x . q {sin 2 [d 2 q ] ■ <S>{z,q) ■ |V(<?)| 2 }(5) 

of the geometric sinusoidal Fresnel term, the fluctuation 
power spectrum, and the visibility of the source surface 
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contains the Fresnel scale from (||). The Fourier transform 
relates the 2-D vector x in the plane of sky, given by the 
vector of the velocity v to the 2-D wave vector q in the 
scattering medium. In (^|) we have neglected a refractive 
cut-off in Fourier space for large q. For quenched scintilla- 
tion the first relevant cut-off in q is in the visibility V(q). 

Anisotropic turbulence can be accounted for in $(z, q) 
but is not considered here. Non-circular and multi- 
component models for the source enter (|^) through an- 
gular variations of the visibility in the q-plane. 

3. The Local ISM 

The ISM of the solar neighbourhood is structured in 
cloudlets, sheets, and walls. The most prominent struc- 
ture is the local bubble (e.g. Breitschwerdt, Freyberg & 
Egger p000| ) with a mean radius of ~ 80 pc. Perpendicular 
to the galactic plane this bubble extends out to about 
200 pc, but its boundary is uncertain due to instabilities 
in the wall surrounding the bubble. The immediate vicin- 
ity of our solar system contains only small amounts of dust 
and a low density ionized phase out to some tens of pc. 
Evidence for absorption in the EUV of white dwarfs gives 
a mean density of 0.04 atoms/cm , an average bubble ra- 
dius of ~ 80 pc, and a five- fold increase i n the gas density 
at the bubble boundary (Warwick et al. 1993| ). 

The amplitude of turbulence in the coronal gas of the 
local bubble is an order of magnitude lower than elsewhere 
in the ISM. This is deduced from scintillation studies of 
pulsars, which show reduced scintillation for the nearby 
pulsar 0950+08, located near the edge of the local bubble, 
relative to other line of sights through the ISM (Phillips 
& Clegg ^992| ). 

While quenched scintillation of extended sources 
favours contributions from the nearby ISM (see ( ]To| ) be- 
low), the scattering material is probably associated with 
the wall of the local bubble. In interactions with neigh- 
bouring bubbles, radio loops are produced and the bubble 
interface is Rayleigh- Taylor unstable. This may have gen- 
erated the large amplitude turbulence in the ionized high 
density gas, necessary for scintillation. For scattering in 
localized regions like the wall of the local bubble the scat- 
tering measure 

rD+H/2 

SM = / AzC 2 N {z) (7) 

JD-H/2 

can be separated in (||) and the medium can be treated as 
a slab of thickness H at distance D. 



4. The Source Structure 

The relation (||) between autocorrelation of intensity fluc- 
tuations and source visibility opens the path to an in - 
tensity interferometer (Hanbury Brown & Twiss 1956 ) 
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Unlike other interferometric techniques in the radio band, 
all phase information is lost, because scintillation in 
quenched scattering does not rely on constructive and 
destructive interference of disturbed wavefronts. The ap- 
pearance of the Fresnel term and the steep power spec- 
trum of turbulence in (||) separates the 'IDV interferom- 
eter' from the optical intensity interferometer and limits 
the angular scale to a small window in 9 at which (^|) can 
be used. 

The time series of intensity fluctuations samples the 
autocorrelation p(x) for different baseline lengths, x, in the 
sky but only along a given direction due to earth motion. 
It is therefore impossible to solve the inverse problem of 
(||) to get the source visibility V(q) from single time series 
of a few days or weeks. 

Even the comparison of source models and the result- 
ing theoretical p(x) with observed IDV is hampered by our 
insufficient knowledge of turbulence in the ISM. Any mea- 
surement of size and structure of the extragalactic source 
is affected by uncertainties in distance, scattering measure 
and relative velocity of the scattering medium. 

From (JSJ) we can derive the square of the modulation 
index m 2 = lim r _»o p(r) of flux variations with a Gaussian 
brightness distribution for the source of width 9: 



p(0) = m 2 



2 \\Def- 2 SM- F(3 ^ /2) (8) 



.DO 2 J " v '' 2 
But it turns out that the structure function 



(9) 



rises for small time-lags as SF oc t 2 independent of the 
spectral index f3. The observed SF show always a much 
shallower rise (see Fig.l) with SM oc r Q ; a = 0.7. . . 1.2 
and are inconsistent with the Gaussian model. The vari- 
ability time-scale of IDV troy derived from the first max- 
imum of SF gives the source size 9 « v iiDv/(l-2 D) inde- 
pendent of the scattering measure. 

Another (extreme) source model is a disk of size 9 with 
constant surface brightness. In this case we find 



p(0) = m 2 = 2 \\D9f- 2 SM ■ F x (/3) 



(10) 



The special functions Fi(/3) & G%{0) are of order unityQ. 
For a succession of slabs at different distances the slab 
closest to the observer dominates the variations according 
to m 2 oc SMDP~ 4 for power spectra (3 < 4 and for sources 
larger than the Fresnel scale in the relevant slab. 



2 The functions Fi (/3) & Gi ((3) axe characteristic for a model 
with constant surface brightness and depend only on the spec- 
tral index f3: 



Fi(/3) = 
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Fig. 2. Autocorrelation of scintillation pattern for a two 
component source model. The components have Gaussian 
surface brightness width 9 — 50/ias. The separation is 79. 



The SF for time-lags much smaller than the correlation 
time iiDV rises in this case as 



SF(r) = 2p(0) ( 
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The characteristic time-scale ^idv °f IDV derived from SF 
gives the source size 

= ^-™ (14) 
D 2 V ' 

in this model. 

The expected intrinsic brightness temperature T^, of 
extragalactic flat-spectrum radio cores is expected to be 
between the equipartition temperature of 4.7 10 10 K and 
the inverse Compton limit of about 6.7 10 11 K. For a given 
Tb we can estimate the distance to the scattering medium 
from ( fl4| ) to be 

-1/2 

(15) 



D = 137pc 
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where we assumed a velocity of v = 20 km/s and trov = 
1 day. 5 is the usual Doppler boosting factor and F v the 
observed flux density of the scintillating component. This 
distance agrees well with the distance to the wall of the 
local bubble. 

So far we have assumed the source to be spherically 
symmetric despite the fact that most IDV sources have 
one-sided jets. 

An indication of deviation from spherical symmetry is 
seen in Fig. 1 where the 'plateau' in the structure func- 
tion is not flat but has at least one minimum beyond 
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Fig. 3. Structure function for a two component source 
model. The source model is like the one in Fig.||. The 
SF are calculated for motion along (solid), perpendicular 
(dashed), and at 14° to the jet axis. 

iiDV- Besides anisotropy of turbulence in the ISM, a two- 
component source structure can lead to this effect as seen 
in Fig. 3. The autocorrelation p{x) of a two-component 
model with two components of 9 = 50/xas each and the 
same flux, separated by O2 — 79 is shown in Fig. 2. p(x) 
has two symmetric secondary maxima along the line con- 
necting the two components in the sky. In such a model the 
minima in SF appear only when the velocity v is within 
±15° of the 'jet-axis' connecting the two components as 
seen from Fig. 3. 

5. Earth Orbit Synthesis? 



Macquart & Jauncey (|2002| ) point out that the degeneracy 
of aligned baselines is broken by changes of earth orbit 
direction. Within 6 months the projected earth motion 
scans 180° of the v ellipse in the sky. In that way the 
plane of p(x) like the one in Fig. 2 can be sampled along 
several directions. In principle the Fourier transform of (^|) 
can then be inverted. Besides knowledge of ISM properties 
as mentioned above, the source must not suffer structural 
changes within at least 3 months. 

In the course of 3 months the distant quasar can un- 
dergo structural changes. Flat spectrum core-dominated 
radio sources are certainly jet-sources and emit boosted 
synchrotron radiation. Structural changes and flux varia- 
tions of the source are limited by the light travel time for 
the source radiu^e^f?. For a source at a medium distance 
z = 1 we get a lower limit on the variability time of 
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which is about 4 months and therefore comparable to the 
required synthesis time. 

3 d,A is the angular diameter distance. For the cosmological 
parameters we use H = 70 km/s Mpc~ x and q — 0.5. 



6. Conclusions 

Quenched scintillation gives a consistent picture of 
Intraday Variability for most IDV sources. The time-scale 
of about 1 day fits with scattering in the ionized gas of 
the wall of the local bubble of a moderately boosted syn- 
chrotron source with brightness temperatures below the 
inverse Compton limit. 

From the analytic expressions for modulation index m 
and structure function at small time-lags presented here, 
it is possible to determine two of the four unknowns: dis- 
tance to the scattering medium, scattering measure, ve- 
locity of the medium, and size of the source. In the slab 
model of the scattering medium together with a constant 
surface brightness of the source it is possible to determine 
the slope (3 ~ 4 of turbulent power spectrum from the 
structure function. 

The structure of IDV sources beyond the size of the 
most compact component cannot be uniquely determined 
from existing experiments. Because most IDV sources 
are intrinsically variable on time-scales of several months 
earth orbit synthesis cannot be reliably realized. 



Extreme scattering events (Fiedler et al. 1987 ) show 
that the ISM is structured on scales of several AU down 



to 0.01 AU (Cimo et al. |200J) well into the scintillation 
domain discussed here. Turbulence in the ISM is there- 
fore not homogeneous on scales of v (3 months) = 1 AU. 
Changes of the long term variability pattern like the sud- 
den end of IDV in 0917+624 (Fuhrmann et al. |2002| ) can 
be due to changes intrinsic to the source or in the scatter- 
ing measure of the ISM. 
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